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ABSTRACT 

Context. Since the discovery, with the EINSTEIN satellite, of strong X-ray emission associated with HD 93162 (=WR 25), this object has 
been predicted to be a colliding-wind binary system. However, radial- velocity variations that would prove the suspected binary nature have yet 
to be found. 

Aims. We spectroscopically monitored this object to investigate its possible variability to address this discordance. 

Methods. We compiled the largest available radial-velocity data set for this star to look for variations that might be due to binary motion. We 
derived radial velocities from spectroscopic data acquired mainly between 1994 and 2006, and searched these radial velocities for periodicities 
using different numerical methods. 

Results. For the first time, periodic radial-velocity variations are detected. Our analysis definitively shows that the Wolf-Rayet star WR 25 is 
an eccentric binary system with a probable period of about 208 days. 

Key words, stars: binaries: spectroscopic - stars: massive stars - stars: Wolf-Rayet - stars: individual: HD 93162 



■ 1. Introduction 
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Massive stars of spectral type O and Wolf-Rayet (WR) are 
important objects that play a crucial role in the dynamic and 
chemical evolution of galaxies. They are the major source of 
ionizing and UV radiation and, through their huge mass-loss 
rates, they have a strong mechanical impact on their surround- 
ings. Despite their importance, our knowledge of these objects 
and of their evolution is still fragmentary. The parameters that 
predominantly determine the evolution of a massive star are 
its mass (as for any star in the HR diagram) and its mass-loss 
rate, although rotation could also have an important impact 
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jMevnet & Maederl [2005). In this context, massive O+O and 
WR+O binaries are key objects because their binary nature al- 
lows us to determine minimum masses from the radial velocity 
orbital solution. 

It is now well established that stron g winds of massive 
stars in binary systems collide (Stevens et al., 1992) gener- 
ating hard X-ray emission. For O stars, this comes in addi- 
tion t o a softer component that is intrinsic to individual ob- 
jects ( Berghof er et all 1 1 997 : San a et all 12006). The existence 
of intrinsic emission for WR stars is still under investigation 
JOskinova et all 120031 iGosset et all [2005). As for colliding- 
wind binaries, the observed characteristics of the collision re- 
gion may contribute to constrain the mass-loss rates and, even- 
tually, the inclination of those systems (by studying the X- 
ray emission modulation during the orbital cycle). Although 
colliding-wind binaries are interesting objects, only a handful 
of them have been studied in detail, which is particularly true 
for WR+O systems. 

The star HD 93162 (WR 25 in Ivan der Huchtl 1200 lb is 
a bright (V = 8.1) galactic Wolf-Rayet located in the Carina 
Nebula region. Its binary nature has been a matter of de- 
bate for man y years. It has been classified as WN7 + 07 by 
smi th ( 1968) because the spectrum appears as a superposi- 
tion of WN-type emission lines and of absorption lines cor- 
responding to an O-type star. The morphology of the blue 
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optic al spectrum of WR 25 has been further discussed by 
Walborn He classified this star as WN6-A (see also 

Walb orn & Fitznatrickl |2000) and considered that there was 



not sufficient information to confirm its binary nature, since 
absorption lines intrinsic to a WN star had been observed in 
the neighbouring binary HD 92740s WR 22 dNiemelaifl973l) . 
The massive binary status of WR 22 with absorption and emis- 
sion li nes moving in phase has been confi r med by several au- 
thors jNie mela. Il973t IConti et all 1 19791: iRauw et all Il996l 
ISchw eickhard t et all Il999|) . In later classification systems of 



Table 1. Observing runs. The first column gives the month and 
year of the observations, the second column lists the instrumen- 
tal combination utilized. The third column gives the number of 
spectra. 



WR stars (e.g. Ivan der Hucht et all fl98ll ISmith et all 1 1996 
the presence of absorption lines has been noted in different 
ways, e.g. "+abs" or adding an "a", indicating that absorp- 
tion lines of unknown origin are observed in the otherwise 
emission-line spectrum. 

Radial velocity studies of WR-type spectra with absorp- 
tion lines are thus needed to shed light on the origin of these 
lines. In the case of WR 25, previously reported radial ve- 
locity studies have not revealed any orbital motion. From the 
study of 15 photog raphic spectra obtained during consecutive 
nights, Moffat ( 1978) concluded that WR 25 is probably a sin- 
gle star. IConti et alJ i 19791) reached a similar conclusion al- 
though they noticed a larger scatter in the radial velocities of 
WR 25 than that observed for other stars in their study (such as 
HD 93 1 3 1 eWR 24 and WR 22 besides its orbital motion). 

Using EINSTEIN observations, ISeward & Chlebowski 

il 19821) found that WR 25 has an abnormally high X-ray flux 
and, later, PolloSi J1991I) suggested that this very high X-ray 
luminosity might be caused by colliding winds in a binary 
system. More recently, iRaassen et alJ J2003h analyzed XMM- 
Newton observations and provided further evidence that the 
X-ray emission could be indicative of a colliding-wind binary 
(CWB) although its apparent stability has been, on th e contrary, 
taken as an argument against binarity. IPollock & Corcoran 
(2006) reported significant variability as detected with XMM- 
Newton and discussed the possibility of periodic variations in 
the X-ray flux. 

Dris sen et alJ \ 19921) found polarimetric variability in 
WR 25, as well as a wavelength dependence of its polarization 
angle. One of the two proposed explanations was that this vari- 
ability may arise in a long-period binary system. The alterna- 
tive explanation attributed particular properties t o the interven- 
ing in terstellar medium. In his recent catalog, Ivan der Huchll 
J200ll) classified WR 25 as WN6h+04f considering that the 
WR emission lines showed evidence of being diluted, although 
no radial-velocity variations had ever been detected. 

To investigate the binary status of WR 25, the massive star 
research groups of Liege and La Plata independently collected 
high-resolution spectra of this star over the past 10 years. The 
acquisition of high-resolution spectra of WR stars is rather un- 
usual due to the difficulties presented by the comparatively 
broad emission lines typical of these stars. However, the radial 
velocities (RVs) measured for each data set showed variations 
larger than the expected errors. No single periodicity could be 
found until the data obtained by the two groups were combined. 
In the following, we show how the analysis of the combined 
spectroscopic data set enabled us to reveal the binary nature 



Date 


Instr. configuration 


n 


High-resolution spectra 


Feb-97 


Echelle-REOSC, 2.15-m, CASLEO 


2 


Mar-97 


CES+LC, CAT 1.4-m, ESO 


4 


Feb-98 


Echelle-REOSC, 2.15-m, CASLEO 


5 


Feb-99 


Echelle-REOSC, 2.15-m, CASLEO 


3 


May-01 


FEROS, 1.5-m, ESO 


4 


Mar-02 


FEROS, 1.5-m, ESO 


4 


Mar-02 


EMMI-REMD/Ech, NTT, ESO 


2 


Apr-02 


FEROS, 1.5-m, ESO 


1 


Jan-03 


Echelle-REOSC, 2.15-m, CASLEO 


4 


Mar-03 


Echelle-REOSC, 2.15-m, CASLEO 


3 


Apr-03 


Echelle-REOSC, 2.15-m, CASLEO 


5 


May-03 


FEROS, 2.2-m, ESO 


3 


Dec-03 


Echelle-REOSC, 2.15-m, CASLEO 


2 


May-04 


FEROS, 2.2-m, ESO 


3 


Apr-05 


Echelle, 2.5-m, LCO 


4 


Feb-06 


Echelle, 2.5-m, LCO 


1 


Medium-resolution spectra 


Feb-94 


2D-Frutti, 1-m, CTIO 


2 


Apr-96 


B&C, 1.5-m, ESO 


2 


Feb-02 


CSpec, 1.5-m, CTIO 


4 


Low-resolution spectra 


May-73 


Cass, 0.9-m, CTIO 


2 


Dec-97 


B&C, 2.15-m, CASLEO 


4 


Apr-01 


REOSC-DS, 2.15-m, CASLEO 


6 



of WR 25 and to derive a preliminary orbital solution for this 
system. 

2. Observations and data reduction 

Our observational data set consists of 50 high-, 8 medium- and 
12 low-resolution spectra obtained, from 1994 to 2006 (except 
for two spectra acquired in 1973), at different observatories and 
with various instrumental configurations. 

Comparison (wavelength calibration) spectra were ob- 
tained at the same telescope positions as the stellar spectra 
immediately after or before the science exposures, except for 
FEROS which is stable enough that three to four calibrations 
over the night are sufficient. 

The summary of the observations is presented in Table ^ 
where we give the dates of the run, the instrumentation used, 
and the number of spectra obtained. 

2.1. High-resolution spectra (HRS) 

Twenty four echelle spectra were obtained with the REOSC 
Cassegrain spectrograph (on long-term loan from the 
University of Liege) attached to the 2.15-m reflector at the 
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Complejo Astronomico El Leoncito (CASLEO 1 ), using a Tek 
1024 x 1024 pixel CCD as the detector. These spectra cover 
an approximate wavelength range from 3600 to 6000 A, with 
a spectral resolving power of 26000. Typical S/N ratios in the 
continuum range between 65 and 100. 

Four spectra were acquired with the ESO 1.4-m Coude 
Auxiliary Telescope feeding the CES spectrograph equipped 
with the Long Camera and CCD ESO#38, yielding an effective 
resolving power of 65000, and a wavelength domain between 
4035 A and 4080 A. Only the blue path was used and the expo- 
sure time was between 30 min and one hour, implying spectra 
with S/N ratios of 100. Normalization to the continuum was 
performed through the use of a metal-poor star observed under 
similar conditions to WR 25. 

Fifteen spectra covering the whole optical range (3750- 
9000 A) were obtained with the Fiber-fed Extended Range 
Optical Spectrograph (FEROS), an echelle spectrograph 
mounted at the ESO 1.5-m telescope at La Silla and then 
transferred to the ESO/MPI 2.2-m telescope in October 2002. 
The detector was a 2kx4k EEV CCD with a pixel size of 
15 /mix 15 jum. The spectral resolving power of FEROS is 
48000. Typical exposure times were between 10 and 20 min 
according to weather conditions, resulting in typical S/N ratios 
in the continuum between 150 and 200 at the 1 .5-m and around 
250 at the 2.2-m. 

Two spectra were acquired with EMMI attached to the 
Nasmyth focus of the NTT telescope. The instrument was used 
in the echelle spectrographic mode (REMD-echelle) with grat- 
ing #9 and grism #3, providing a resolving power of 7700 
and 18 usable spectral orders, covering the wavelength domain 
from 4040 A to 7670 A. Typical S/N ratios are around 120. 

Five spectra were observed with the 2. 5-m du Pont tele- 
scope at Las Campanas Observatory (LCO), Chile, using the 
echelle spectrograph which provides simultaneous wavelength 
coverage from ~ 3700 to 7000 A at a typical resolving power 
of -40000. The detector is a Tek5 2k x 2k CCD with pixel size 
of 24 /mi. We used a 1 x 4 arcsec slit. These spectra have S/N 
ratios around 100. 

2.2. Medium-resolution spectra (MRS) 

Two medium-resolution spectra (S/N ~ 100) were obtained in 
1996 with the ESO 1.5-m telescope equipped with the mod- 
ified Boiler & Chivens spectrogr aph at La Si l la. Th e con- 
figuration utilized is described in iGosset et al.l J200lh . Four 
medium-resolution spectra were obtained with the Cassegrain 
Spectrograph at the CTIO 2 1.5-m telescope in February 2002. 
A 600 1/mm grating blazed at 3375 A in second order, com- 
bined with a Loral lk CCD provided a 2 pixel resolution of 
1.5 A (wavelength range 3270-4180 A). Two spectra were ob- 
served with the Shectman/Heathcote two-dimensional, photon- 
counting detector (2D-frutti) on the Cassegrain spectrograph at 

1 CASLEO is operated under agreement between CONICET, 
SECYT, and the National Universities of La Plata, Cordoba and San 
Juan, Argentina. 

2 CTIO is operated by the AURA Inc. under a cooperative agree- 
ment with the National Science Foundation as part of the NOAO. 



the CTIO 1-m telescope. The wavelength coverage was from 
3800 to 5000 A at a 3 pixel resolution of 1.5 A. Typical S/N 
ratios range between 50 (2D-frutti data) and 250 (CSpec at 
CTIO). Typical resolving powers for these MRS are in the 
range 2000-4000. 

2.3. Low-resolution spectra (LRS) 

We have also obtained lower resolution spectra with resolving 
powers ~ 1000. 

Six LRS were acquired at CASLEO with the above- 
mentioned REOSC spectrograph but in its simple dispersion 
mode. This configuration provides a sampling of 1.64 A per 
pixel, on the Tek 1 024x 1 024 CCD. 

We also used 4 spectra observed with the Boiler & Chivens 
spectrograph attached to the 2.15-m telescope at CASLEO. A 
PM 512x512 pixel CCD, with pixel size of 20 /mi, was used as 
the detector. The reciprocal dispersion was ~ 2.3 ApixeL 1 , and 
the wavelength region observed was about AA 3800 - 4800 A. 
Typical S/N ratios of CASLEO LRS are about 200. 

In addition, we used two previously unpublished radial ve- 
locity values determined from spectra observed in May 1973 at 
CTIO with the Cassegrain spectrograph attached to the 0.9-m 
telescope. These spectra were recorded on photographic plates, 
have a reciprocal dispersion of 120 A/mm, and were widened 
to 1 mm for a better visibility of the spectral lines. The spec- 
trograms were measured for the determination of RVs with a 
Grant oscilloscope engine. Because the radial velocity values 
of the Niv /14058 emission determined from these two spectra 
were rather more negative than the mean of the other 52 ob- 
servations, they were deemed suspect and were not included in 
the radial velocity study published by IConti et alJ Jl979l) . We 
note, however, that the radial velocity of the interstellar Ca n 
absorptions in the two photographic spectra (-32 km s~') does 
not deviate much from the values for these lines in our digital 
spectra (see below). 

3. The radial velocity analysis 

We have determined the radial velocities of WR 25 measur- 
ing the position of the N iv /14058 emission line because this 
line is narrow and strong enough to minimize measurement er- 
rors. In addition, it is expected to be formed relatively deep 
in the WR wind, and thus to better reflect the true motion 
of the WR star. We arbitrarily adopt ed for that l ine th e rest 
wavelength /t =4057.9 A as given bv lConti etail Jl977l) . The 
adopted method consists of fitting Gaussian profiles to the ob- 
served line using either MIDAS or IRAF routines. We worked 
on spectra normalized to the continuum and, to homogenize the 
RV measurements, we favoured the position as measured on the 
upper part of the line; this has the advantage of being less de- 
pendent on the errors in the definition of the continuum. The 
RVs of the N iv /14058 emission line as measured in the spectra 
of WR 25 are reported in Table[2] These values should be cau- 
tiously considered as preliminary because various effects (e.g. 
intrinsic line-profile variations) could render the measured RVs 
inaccurate. On a few good spectra, we used other techniques to 
measure the position of the Niv /J4058 line, e.g. one based on 
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Table 2. Radial velocities (RVs) of the Nrv ,14058 emis- 
sion line measured in the spectra of WR 25. The velocities 
are expressed in the heliocentric rest frame; time is given in 
Heliocentric Julian Date. 



HJD 


RV 


Dataset 


HJD 


RV 


Dataset 


2,400,000+ 


km s" 1 




2,400,000+ 


km s _I 




41823.657 


-112 


LRS 


52328.707 


-18 


MRS 


41828.620 


-117 


LRS 


52328.718 


-20 


MRS 


49406.668 


-9 


MRS 


52333.848 


-20 


MRS 


49409.767 


-23 


MRS 


52335.652 


-18 


HRS 


50181.671 


-32 


MRS 


52337.685 


-27 


HRS 


50182.637 


-25 


MRS 


52338.636 


-20 


HRS 


50505.846 


-47 


HRS 


52339.639 


-22 


HRS 


50507.802 


-48 


HRS 


52353.616 


-20 


HRS 


50531.535 


-68 


HRS 


52353.620 


-20 


HRS 


50531.584 


-70 


HRS 


52383.507 


-48 


HRS 


50532.513 


-59 


HRS 


52655.865 


-29 


HRS 


50534.524 


-69 


HRS 


52657.870 


-29 


HRS 


50807.849 


-15 


LRS 


52658.861 


-30 


HRS 


50809.868 


-17 


LRS 


52659.855 


-29 


HRS 


50811.869 


-28 


LRS 


52710.711 


-4 


HRS 


50812.856 


-26 


LRS 


52711.827 


-6 


HRS 


50842.800 


-7 


HRS 


52712.832 


-9 


HRS 


50846.791 


-7 


HRS 


52735.672 


-11 


HRS 


50847.877 


-7 


HRS 


52736.673 


-18 


HRS 


50850.855 


-11 


HRS 


52737.631 


-15 


HRS 


50852.798 


-8 


HRS 


52738.590 


-12 


HRS 


51209.865 


-15 


HRS 


52739.658 


-22 


HRS 


51216.881 


-5 


HRS 


52782.493 


-32 


HRS 


51218.880 


-10 


HRS 


52783.498 


-30 


HRS 


52005.575 


-80 


LRS 


52784.493 


-30 


HRS 


52007.577 


-102 


LRS 


52985.841 


-42 


HRS 


52008.576 


-110 


LRS 


52989.842 


-37 


HRS 


52009.612 


-81 


LRS 


53131.510 


-11 


HRS 


52011.635 


-103 


LRS 


53133.605 


-10 


HRS 


52013.589 


-74 


LRS 


53135.485 


-11 


HRS 


52037.638 


-16 


HRS 


53480.710 


-46 


HRS 


52038.586 


-20 


HRS 


53481.584 


-39 


HRS 


52039.618 


-22 


HRS 


53488.559 


-28 


HRS 


52040.622 


-15 


HRS 


53490.662 


-26 


HRS 


52328.697 


-22 


MRS 


53772.679 


-12 


HRS 



the lower part of the emission profile. The resulting RVs differ 
from the adopted RVs measured on the basis of the upper part 
by 0.06 A or 4.5 kms -1 at maximum. No systematic shift has 
been detected. These results give further support to the reported 
RVs and provide an idea of the error on the measurement. 

In Fig. ^ we show two typical line profiles. One was ob- 
served on HJD 2452 383.507 (corresponding to cp = 0.78, see 
below) and exhibits a RV of -48 km s _1 . The other, acquired on 
HJD2453 133.605 (<p = 0.39), corresponds to a velocity of -10 
kms -1 . These spectra are representative of two extreme posi- 
tions of the line as observed by us at high resolution. The main 
difference between the two spectra of Fig.^is a Doppler shift 
and line-profile variations, if they exist, are only second order 
effects. On the other hand, the smallness of the shift compared 
to the line width explains the difficulty in performing fully ac- 
curate measurements through a fit of the entire profile. 



a L6 

o 

m 
-a 

N 1.2 



1.4 



o 



1 1 1 1 1 1 II 


1 1 1 1 1 1 1 1 1 1 


. 1 1 1 1 1 

= 0.78 _ 






= 0.39 . 








I I I I 1 I I I I 1 I I I I 1 I I I I 1 I I I 



4045 4050 4055 4060 4065 4070 

Wavelength(A) 

Fig. 1. Two typical line profiles corresponding to high S/N 
FEROS spectra of WR 25. The emission line is N iv /14058 and 
is shown in the heliocentric rest frame. One (plain line) was ob- 
served on HJD 2452 383.507 (corresponding to = 0.78, see 
Table [3j and is exhibiting a RV of -48 kms -1 ; the other was 
observed on HJD 2 453 133.605 (0 = 0.39) corresponding to a 
velocity of -10 km s _1 . The main change between the two spec- 
tra is a shift in RV. 



We also measured the Can /13933 interstellar absorption 
line. In the echelle spectra, this line presents at least 5 distinct 
components. Therefore, we measured the position of the central 
line of the three main components. Averaging the RVs obtained 
in the echelle spectra, we derived a mean of -30±4kms _1 . In 
Fig- El (bottom), we also show the RVs of this Can line mea- 
sured in our WR 25 spectra. This confirms the good agreement 
existing among the different instrumental configurations used 
in this work. 

Having noticed that the position of the N iv /14058 line was 
variable, we studied the time series of the measured RVs. We 
searched for periodicities using two independent methods: the 
algorithm to derive period s of cyclic phenomena described in 
Marraco & M uzzio](ll980l) . and a Fourier-type analysis method 
by H eck et alJ(ll985l see also comments by Gosset et al., 2001). 

We first analyzed the data set consisting of the 50 high- 
resolution spectra (HRS). The particular distribution of the ob- 
serving times induces some aliasing at 5v = 0.00045 d , at 
5v = 0.00270 d -1 , at 5v = 0.02540 d 1 and at 6v = 0.02810 d _1 
as revealed by e .g. the spectral window (see the definition in 
lDeemingi fl975]). The Amplitude Spectrum (square root of the 
power spectrum) of the HRS data is given in Fig. 13 The highest 
ordinate is located at v = 0.004795 d _1 . The two neighbouring 
peaks are aliases of this frequency; this is also the case for the 
peaks at v = 0.00209 d _1 , at v = 0.00749 d _1 , at v = 0.03015 d" 1 
and at v = 0.03286 d _I . All of them thus belong to the same 
family and we consider that the dominant frequency is the one 
at v = 0.004795 d _I . This is further confirmed by a decomposi- 
tion of the RV curve into several frequencies u sing the mul- 
tifrequency approach (see Eqs. All to A19 in Gosset et"aQ 



2001). The main frequency varies from v = 0.004795 d" 1 to 
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Fig. 2. Amplitude Spectrum (square root of the power spec- 
trum) as a function of the frequency expressed in d -1 and cor- 
responding to the analysis of the HRS data set. From 0.04 d" 1 
to 0.5 d _I , no peak exceeds 13 kms~'. The adopted main fre- 
quency is marked by a filled circle and the expected positions 
of the main aliases of this frequency are marked with open cir- 
cles. 



v = 0.004805 d 1 depending on the secondary frequencies be- 
ing inventoried: we thus adopt v = 0.00480 d" 1 (P = 208.3 d) 
as the true progenitor. We estimated the error on the period 
by taking as a conservative upper limit one tenth of the peak 
width (which is a function of the time baseline); we derived 
a-y = 0.00003 d" 1 and thus <7> = 1.3 d. The Marraco & Muzzio 
method applied to the RVs measured on the HRS spectra gave 
as a most probable period P = 208.3+0.5 d along with some 
aliases, in good agreement with the Fourier results (see Fig.|3}. 

A phase diagram with the adopted period presents a clear 
gap (around phase 1.0, see Fig0}. Therefore, we complemented 
our HRS data with other data sets. We thus defined two subse- 
quent data sets: 

- our 50 RVs measured on the high-resolution spectra (HRS) 
complemented by the 8 medium-resolution data (MRS) and 
by the 12 low -resolution spectra (LRS), thus 70 data points 
labelled HRS+MRS+LRS, 

- all our 70 RVs plu s the 54 published bv lMoffaddl978h and 
IConti et all dl979h - comple mented with the RV m easured 
in one spectrum provided by Ham ann et al.l 11995). 

We also ran the period-search algorithms using these 
enlarged data sets. In the HRS+MRS+LRS data set, the 
Fourier analysis favoured a frequency v = 0.004820 d _1 
which is within the error bars of the previously quoted 
one. Finally, the data set all+published indicates a frequency 
v = 0.004810 d- 1 (cr v = 7.7 I P" 6 d" 1 ) corresponding to 
P = 207.9 d (cr P = 0.3 d). The IMarraco & Mnzziol |l980) 
method applied to the all+published data set similarly gives 
P = 207.8+0.3 d. 
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Fig. 3. Run of the IMarraco & Muzziol (^980) statistic as a 
function of the frequency expressed in d -1 and correspond- 
ing to the analysis of the HRS data set. The main dips are 
seen at v = 0.0012 d~\ v = 0.0016 d _1 , v = 0.0024 d" 1 
and v = 0.0048 d _1 ; the first three are subharmonics of 
v = 0.0048 d _I and are a well-known artifact of this 
kind of method. Other dips are visible at v = 0.0082 d -1 , 
v = 0.0100 d"\ v = 0.0302 d~ ! and v = 0.0329 d _1 . The last 
two are also aliases of the v = 0.0048 d -1 dip. The cited fre- 
quencies are marked as in Figure|5] 



To further refine the period determination, we ran GBART, 
an improved v ersion 3 of the orbital solu tion program originally 
published by Bertiau & Gro bbenldl968l) . with the three above- 
mentioned data sets. The RVs measured on HRS were weighted 
with 1, MRS with 0.5, and LRS or published, with 0.1. The or- 
bital parameters obtained for each data set are explicitly given 
in Tableland the corresponding orbital solutions are depicted 
in Fig.0] We obtained quite eccentric orbits with e ranging be- 
tween 0.35 and 0.5, orbital semi-amplitudes K between 33 and 
44 km s _1 , and different values for the orbital period. 

The fact that there is no HRS RV more negative than 
-70 km s _1 , while some LRS RVs indicate -1 10 km s _1 , is in- 
triguing. Clearly the HRS data set presents a strong gap around 
phase 1 .0. This gap can be filled by LRS data obtained in May 
1973 and April 2001. Even though the errors on these RVs are 
expected to be larger and despite the low weight given, these 
LRS data have a large impact on the eccentricity and on the 
orbital semi-amplitude derived. Adding the published data to 
ours still yields an orbital solution that remains in good agree- 
ment with the HRS+MRS+LRS one. However, none of the 
published RVs have been acquired near the critical phase (1.0). 

Although the exact orbital solution remains uncertain, our 
various data sets allow us to conclude that WR 25 presents RV 
variations with a period of about 208 days. These variations are 
indicative of an eccentric binary system. Strictly periodic RV 
variations with such a long time-scale are difficult to interpret 



3 Available upon request from 
ftp : //lilen . f caglp . unlp . edu . ar 
/pub/fede/gbart-0. 1-41. tar. gz 
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Table 3. Orbital solutions corresponding to the RVs of the N iv 
/14058 emission line within different data sets. The quoted er- 
rors correspond to lcr uncertainties. Symbols have the canon- 
ical meaning. The last three correspond, respectively, to the 
mass function, the standard deviation of the fit and the num- 
ber of data points involved. 

HRSonly HRS+MRS+LRS All + published 



P[d] 


208.3 


± 0.2 


207.9 


+ 


0.1 


207.85 


+ 


0.02 


V [kms- 1 ] 


-30.6 


± 0.7 


-34.2 


+ 


0.7 


-34.6 


+ 


0.5 


K [kms" 1 ] 


33 


± 2 


42 


+ 


2 


44 


+ 


2 


e 


0.35 


± 0.03 


0.48 


+ 


0.02 


0.50 


+ 


0.02 


at [degrees] 


227 


± 4 


216 


± 


3 


215 


+ 


3 


Tperiast [4]* 


1598 


± 3 


1597 


± 


2 


1598 


+ 


1 


TRVmax [4]* 


1654 


± 3 


1655 


+ 


2 


1654 


+ 


1 


a sin; [Rq] 


125 


± 8 


151 


+ 


9 


156 


+ 


8 


F(M) [Mo] 


0.6 


± 0.1 


1.1 


+ 


0.2 


1.2 


+ 


0.2 


<t [km s" 1 ] 


2.5 




3.7 






4.7 






n 


50 




70 






124 







Heliocentric Julian Date 2,450,000+ 



in other terms. The exact values for some of the parameters are 
highly dependent on low resolution data. Therefore, we need 
further good quality high-resolution data in order to derive a 
fully definitive orbital solution for the WR primary. 

We also considered other lines in the HRS spectra of 
WR 25. Clearly, certain lines such as N v 114603-4619 vary in 
agreement with the derived orbital solutions. However, some 
absorption lines such as Hei ,15876 and ,14471 seem to vary in 
anti-phase, plausibly revealing the signature of the binary com- 
panion. 

Studying other lines will help to verify whether the RVs de- 
rived from the motion of Niv 14058 represent the orbital mo- 
tion of the primary; this can not be established on the basis of 
a single line. This is beyond the scope of the present paper and 
requires high-resolution high S/N data covering all the phases. 



4. Conclusions 

We have found that the RVs of the N iv /14058 emission line in 
the spectrum of WR 25 show variations with a period of about 
P=208 days. These variations indicate that the WN star WR 25 
is an eccentric binary system. 

The He i /15876 and -14471 absorption lines seem to show 
an anti-phase motion and should thus partly belong to the com- 
panion which is supposed to be of the OB type. 

The companion seems to be fainter than the WR primary, 
indicating that WR 25 is perhaps another example of a young 
binary system where the WR component is the most massive 
star. The prototype o f this class is WR 22 (WN7h+09III-V, 
van derHuchdEooi which turned out to be a massive binary 
system (55 M n + 21 M n ) in an 80- day period orbit l lRauw et all 
1 19961 ISchweickhardt et all Il999h . The primary of WR 22 is 
most probably an example of the long-searched for case of a 
core H burning star resembling a WR star due to its luminosity 
and mass. The star is on its way to becoming a hydrogen-free 




Orbital Phase 

Fig. 4. Observed radial velocities (along with the fitted mod- 
els) corresponding to the Niv 14058 emission line in the spec- 
tra of WR 25. The data are phased according to the ephemeris 
as given in Table |5J phase 0.0 corresponds to periastron. The 
different panels illustrate different analyzed data sets. Top: Our 
HRS data only (filled circles). Middle top: HRS+MRS+LRS 
data; open triangles represent the RVs in the MRS set, and 
squares the LRS ones. Mid dle bottom : All + publish e d data ; 
cro sses are the RVs from IConti et all ril979h . Dvloffad Jl978t> . 
and lHamann et all (119951) . In each plot, the solid curve repre- 
sents the relevant orbital solution as given in Table [3] Bottom: 
Radial velocities of the interstellar Can K absorption feature 
measured on our spectra. 



WR-type star. These objects are extremely rare and thus inter- 
esting to study in detail. 

WR 25 is one of the brightest WN stars in the X-ray do- 
main, too bright to be explained by intrinsic X-ray emission. 
It thus was customarily classified as a putative colliding-wind 
binary. Fig.|4]is the first proof that WR 25 is a binary as sus- 
pected from its X- ray emission. On the basis of the existing 
XMM -Newton data. IPollock & Corcoran! (Eo06) found that the 
X-ray emission of WR 25 is variable although the period can- 
not be derived independently from the X-ray data alone. WR 25 
was brighter in the X-ray domain on JD 2452 842.6 which ac- 
cording to our ephemeris corresponds to a phase of 0.96-0.98. 
The star thus seems to be brighter at periastron as predicted 
for colliding- wind binaries in the adiabatic regime (see e.g. 
IStevens et all Il992l) . Although this should be further investi- 
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gated, WR 25 appears to be an example of a colliding-wind 
binary system. 

We are currently planning to monitor the optical spectrum 
of WR 25 around the expected time of minimum radial veloc- 
ity of the WR emission lines in order to derive a full definitive 
orbital solution. High-resolution high S/N ratio data are essen- 
tial. These spectra will help to improve the orbital solution, 
and to detect the absorption lines of the companion, which for 
this particular phase domain should be located to the red of the 
WN emission line, thus avoiding confusion with the absorption 
components of P-Cygni profiles from the primary. 
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